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Abstract 

An atmospheric general circulation model (AGCM) with idealized and complete physics has 
been used to evaluate the Tropical Easterly Jet (TEJ) jet. In idealized physics, the role of upper 
tropospheric friction has been found to be important in getting realistic upper tropospheric zonal 
wind patterns in response to heating. In idealized physics, the location and strength of the TEJ 
as a response to Gill heating has been studied. Though the Gill model is considered to be widely 
successful in capturing the lower tropospheric response, it is found to be inadequate in explaining 
the location and strength of the upper level TEJ. Heating from the Gill model and realistic upper 
tropospheric friction does not lead to the formation of a TEJ. 
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AGCM 


1 Introduction 


The Tropical Easterly Jet (TEJ) is an importan t uppe r tropospheric phenomena that occurs during the 
Indian Summer Monsoon Season (jKoteswaraml (jl958l ). iFlofml (J1968I )). It has been shown that its loca¬ 


tion, s tructu re and strength depend on the location and magnitude of heat sources ([Rao and Srinivasan 


1 2013 . 


20151) ). But the methodology of using an atmospheric general circulat ion model ( AGCM) has 


its own complications because of the detailed physics and parametrizations. iGilll (|198 (T) . henceforth 


G80, showed in his classic paper that heat sources in the tropics directly influence lower tropospheric 
wind and that the patterns bore qualitative resemblances to observations. It is from his paper that we 
know that a jet-like structure is possible even in a simple setting. 

One of the important drawbacks of this model is that it is essentially a linear single-layer shallow 
water model, but actually represents a two layer baroclinic system. The upper and lower level winds 
are viewed as a single response. The zonal wind, which was a response to diabatic heating, was a cosine 
function of height which caused it to reverse signs in the lower and upper troposphere with maxima at 
top and bottom. Both upper level heating and boundary layer convergence are thus tightly coupled. 
It thus ignores the rol e of latent hea ting in driving the surface winds, which is true only if heating is 
at much higher levels (INecli n (|1989l) ). 

The requirement of a strong equivalent linear mechanical damping to get realistic wind response 
is also a consistent critique of this model. The frictional parametrization is identical i n the lower and 
upper troposphere. The source of this damping in the upper troposphere is not clear. 


Lin et al 


(2008!) 


showed that the strong damping used in Matsuno-Gill-type models was spatially inho mogeneous, and 


depen ded on the heating rate, as prescribed by the deep convective or shallow clouds. iGregorv ct al 


(1997) found that the Tropical Easterly Jet in the Met Office (UKMO ) GCM was weakened due to 


damping caused by Convective Momentum Transport parametrization. 


Showma n and Po lvan i ( 2010!) 


showed that in Matsuno-Gill models e quatoria l easterlies pe r sist d ue to the presence of momentum 
transport by the imposed heat source. iHoskins and R.odwelll (11995 1 used a primitive equation model 
that used Rayleigh damping which increased as one moved into the stra tosphere and beyond . The 


use of a multi-layered atmosphere model for studying the TEJ is not new. 


Mishra and Tando nl ( 1983! ) 


showed the importance of beta effect for wave growth of the TEJ. The effects of non-linearity on the 


TEJ was also studied. 


Mishra and Tandonl (119831) : iMishral (11987!) further showed that non-linearity 


was important in making the jet have a longer zonal extent. 

In this paper, an AGCM, the Community Atmosphere Model, version-3.1 (CAM-3.1) has been 


used to study the effect of upper tropospheric friction on the TEJ. Heating as specified from the 


Gill 
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(1980) model has been applied to study the location and structure of the TEJ. The physics in CAM- 


3.1 has been considerably simplified by totally removing the entire default physic s parametrizations. 
Instead the model has been run with specifications as in IHeld and Suarezl (119941 ) (hereafter referred 
to as HS94), unless explicitly stated. Such simulations with idealized physics have been termed as 
‘ideal-physics’. 


2 Model description and experimental details 


The AGCM that has been used is CAM-3.1. The finite-volume (FV) dynamical core with 2°x2.5° 
horizontal resolution and 26 vertical levels, has been used. The time step is the default value of 30 
minutes. The physics is as in HS94, unless explicitly stated. Since CAM-3.1 has hybrid-sigma co¬ 
ordinates, frictional damping from levels 1 to ~0.7 is identical to HS94. Newtonian relaxation to a 
specified equilibrium temperature is also identical to HS94. 

The effect of upper tropospheric friction has been explored by comparing simulations with and 
without upper tropospheric friction, against simulations with full AGCM physics, henceforth called 
full-physics simulations. Simulations had been conducted with the entire physics of CAM-3.1. The 
actual heating (inclusive of both convective and radiative) and temperatures from these simulations 
have been prescribed in the ideal-physics mode of CAM-3.1. Then the impact of upper tropospheric 
friction has been assessed by comparing the location and strength of the jet in ideal-physics simulations 
against simulations with full AGCM physics. 

The full-physics simulations considered are: 


1. Aqua-planet simulation with a single heat source. Precipitation was induced by elevating tem¬ 
peratures with a peak value at 90°E,20°N. The peak value was 32°C which was brought down to 
a value of 25°C. The elevated SSTs had an oval shape with a 90° major-axis as 16° minor-axis. 
This simulation, named AP_90c20n, has been run for twelve months and the data of the 12 th 
month has been used. 


2. A one year simulation of the default configuration of CAM-3.1 with clim atological SSTs, but 


witho ut orography. Since the TEJ is the strongest in the month of July (jRao and Srini vasa n 


(]2013l )), the month of July has been considered. This simulation has been named noGlOrog. 


The names of the corresponding ideal-physics simulations has ‘_ip’ appended to the original names. 
The inclusion of upper tropospheric friction is understood by further appending ‘_nf’ (no friction) or 
‘_pf’ (friction present). For example, for the aqua-planet simulation with upper tropospheric friction, 
the name for the corresponding ideal-physics simulation is AP_90e20n_ip_pf. The temperatures from 
the full-physics simulation have been zonally averaged and then imposed in the ideal-physics cases. 
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Rayleigh friction is imposed in the upper tropospheric friction. Damping, which is a function of height, 
was initiated at a hybrid level of of 0.3. A 2 nd degree polynomial was chosen to fit values at three 
data points: (i) damping timescale of 30 days at 0.25 hybrid level, (ii) damping timescale of 300 days 
at 0.085 hybrid level and (iii) damping timescale of 1000 days at 0.0 level. Table [T| gives details about 
the simulations. In the 3 rd column, CAM-3.1 implies that frictional parametrization is the default as 
in full-physics of CAM-3.1. The ideal-physics simulations have been run for 12 months and the last 
six months have been averaged to study the results. 

3 Effect of upper tropospheric friction 

Before proceeding, a brief description of some terms are given. Q max denotes the maximum heating 
rate. U ma x denotes the value of peak zonal easterly and U norm denotes the zonal velocities normalized 
by U m ax- In the figures U norm contours are shown to better facilitate comparison between full-physics 
and ideal-physics simulations. Since the winds are normalized by U maxi the positive contours in the 
figures are actually easterlies, and the negative contours are westerlies. Jet length and width are 
defined on basis of the location of U max and closure of the U norm = 0.5 contour, that is, closure of 
the 0-5U max contour. Arcs in zonal and meridional directions passing through the location of U max 
are considered. Zonally and meridionally, the length in degrees when these arcs touch the 0.5 U max 
contour are considered as jet length and widths respectively. If the 0.5 U max contour does not close 
zonally, then the jet is defined to be non-existent. 

3.1 Aqua-planet simulation 

The aqua-planet configuration used has been described in section[l]above. The maximum heating rates 
is shown in Fig. Hal The maximum heating has been shown because this enables easy understanding 
of peak heating, since at different locations the rates peak at different pressure levels . Table Q] has 
details about the full-physics and corresponding ideal-physics simulations. 

For comparison purposes the U norm contours of actual aqua-planet simulation of the 12 th month 
is shown in Fig. Dal For the ideal-physics simulations, the differences were found to be minor when 
either temperatures at each grid point or zonally averaged temperatures were specified. Hence we only 
show the latter. This lack of difference demonstrates that provided radiative-convective equilibrium 
temperatures are prescribed, heating plays a dominant role in determining the profile of the TEJ. 

In the absence of friction, it is seen from Fig. [2c] that the zonal location of U max is eastwards by 
20° compared to full-physics simulation (l2al) . However, use of upper tropospheric friction causes the 
location of U max (Fig. O to be very similar to full-physics simulation. Table 0 confirms this. Similar 
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(a) AP_90e20n (SST: 32°C peak on 25°C background) (b) noGlOrog (no orography) 


Figure 1: Maximum heating ( K day 1 ) rate of AP_90e20n and noGlOrog simulation 


Table 1: Magnitude (ms 1 ) and location of peak zonal wind ( U max ) 


Case 

Upper tropospheri 

Simulation type 

friction 

LC 

Zonal wind 


Peak 

Lon 

Lat 

Press (hPa) 

AP_90e20n 

full-physics 

CAM-3.1 

41.64 

47.5°E 

4°N 

150 

AP_90e20njp_nf 

ideal-physics 

Absent 

42.35 

27.5°E 

0° 

150 

AP_90e20n_ip_pf 

ideal-physics 

Present 

39.36 

47.5°E 

2°N 

150 

noGlOrog 

full-physics 

CAM-3.1 

42.34 

30.0°E 

10°N 

125 

noG10rog_ip_nf 

ideal-physics 

Absent 

46.09 

55.0°E 

6°N 

125 

noG10rog_ip_pf 

ideal-physics 

Present 

44.09 

55.0°E 

8°N 

125 

AP_90e20n: aqua-planet; noGlOrog: no orography 
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trends were also observed in the meridional profiles of U norm . The zonal profile with friction was also 
more realistic compared to the simulation without friction. 

Though not shown, ideal-physics simulations were also conducted using zonally averaged tem¬ 
peratures from noGlOrog simulation described in section [Q In the absence of friction the jet in 
ideal-physics became more zonally extensive (Fig. Eel) . But presence of upper tropospheric friction 
made the zonal profile similar to Fig. [2e]. This can be understood from the thermal wind equation 
(write it). The zonal thermal wind (ut) is proportional to meridional temperature gradient as follows: 
u T = Su g ~ jp (5P), where, R a is gas constant for dry air, / is Coriolis parameter, P is 

pressure and T is temperature. Compared to aqua-planet, heating in noGlOrog is elevated over a 
significant zonal stretch (refer Fig. SB- Thus using zonally averaged temperatures from noGlOrog 
implies that temperatures are on an average elevated in the zonal direction in comparison to aqua- 
planet simulation. Thus meridional temperature gradients do not change over a significant zonal scale. 
This implies a tendency for the zonal wind to be uniform over a greater zonal span. In the absence 
of upper tropospheric friction this leads to a significant increase in jet length. With the location of 
U m ax similar to Fig. [2c] But friction again caused the magnitude of U max and zonal wind pattern to 
be similar to Fig. [Ml This further shows the importance of upper tropospheric friction. 

3.2 noGlOrog simulation 

The heating imposed in the ideal-physics simulations is from July 1 st year full-physics simulation. The 
heating is shown in Fig. Ilbl Fig. [2b] shows the horizontal cross-section of full-physics simulation. As 
with the ideal-physics aqua-planet simulations, ideal-physics simulations were conducted with zonally 
averaged temperature profiles and temperature specified at each grid point. Since was no major 
difference in the structure of the TEJ between the two sets, only the former cases are shown. 

From Figs. [ME we can see that ideal-physics simulation has correctly reproduced the horizontal 
structure of the TEJ. Although the location of U ma x is much westwards compared to full-physics, 
there is a broad region of zonal wind maxima. The velocity differences in the region enclosed by the 
0.9 contour was found to be about 2-3 m s -1 . This region is more zonally extensive than the full- 
physics simulation. It is also seen that upper tropospheric friction does not have much influence on 
the simulation except that it slightly reduces the value of U ma x (refer Table [T]). This further implies 
that when, without upper tropospheric friction, the ideal-physics simulation is realistic, presence of 
friction does not make it unrealistic. 
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(a) AP_90e20n: full-physics (b) noGlOrog: full-physics 



(c) AP_90e20n_ip_nf: ideal-physics 



(d) noG!Orog_ip_nf: ideal-physics 



60°W 30°W 0° 30°E 60°E 90°E 120°E 150°E 180° 

(e) AP_90e20n_ip_pf: ideal-physics 



(f) noG!Orog_ip_pf: ideal-physics 


Figure 2: U norm contours. Left panel is AP_90c20n simulation set, right panel is noGlOrog simulation 
set. The ideal-physics simulations have heating and zonally averaged equilibrium temperatures from 
the corresponding full-physics simulation. Cross-sections are at U ma x pressure level, ‘cross-diamond’ 
is location of U ma x- Grey region denotes easterlies. Additional details are in Table [l] and Fig. [I] 
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4 The Tropical Easterly Jet in Gill model and its comparison 


with CAM-3.1 


4.1 Zonal wind from Gill model 


There has been a lack of interest in study ing the cha racteristics of the Tropical Easterly Jet. This 


Gill (I1980f) which is a linear shallow water mode l using 


is partly due to the remarkable paper by 
the equatorial /3-plane approximation. Although critical analysis and criticisms (for example 
(1198911 . Plumb (|2010r il for this model have been given, it still remains a pioneering work. 


Neelin 


The asymmetric solution of Gill’s formulation has been considered. The analytical expressions for 
heating and the resultant zonal winds are given by Eqs. m and ©• 


Q(x,y) = F(x)y exp - 


y 


(1) 


where, 


where, 


F{x) = 


cos (kx) if |x| < L 
0 if \x\ > L 


k = — 

2 L 


93(e,x) = < 


1 / y 2 
i{x,y) = -q 3 (e,x)(y 3 - Gy) exp — 


—k [l + exp (—lOeL) ] exp {5e(a; + L)} if x < —L 

— 5e cos (kx) + k [ sin (kx) — exp (5e(x — L)} ] if |x| < L 
0 if x > L 


e = 0.1, a factor used to represent Rayleigh friction and Newtonian cooling, 

L: non-dimensional zonal length of the domain across which forcing is considered 


( 2 ) 


Qmax = +V2e 2 at (x,y) = (0, y/2) 
Qmin = -V2e~i at (x,y) = (0, -V2) 


The heating and zonal wind reach their extremum values at the locations given below: 














U min = -2.45 at (x,y) = ^-1.14, (6 - 2\/6) 2 ^ 
u max = +2.45 at (x,y) = (-1.14, - (6 - 2^6) 

The three-dimensional structure of heating is obtained by multiplying Eqs. © and m as follows: 


( 7TZ \ 

u) 

u(x, y, z) = u{x, y) cos 


(3a) 

(3b) 


where, 0 < z < D, D being the non-dimensional height of the domain 

Using the standard formulae given in G80, the equations have been non-dimensionalized as follows: 


c = ^fgH, L d 



^ 6d T d 


where, g = 9.81 m s -2 , H = 400 m, c = 62.64 m s -1 (speed of gravity waves in the absence of rotation), 
L d = 10° (horizontal length scale), r d = 4.93 hour (time scale), e d = 0.021 hour -1 (the dimensionalized 
dissipation rate). This then gives the heating to be dimensionalized by Qd=r c j -1 = 4.87 K day -1 . 

Additional details are given in Table[2] The heating maxima is centered at 0°,10°N. The meridional 
and horizontal structure of the normalized zonal winds is shown in Figs. I3al and l3bl respectively. The 
default value of e, 0.1, gives very high zonal winds speeds. For this case the zonal distance of U ma x is 
~11.5° to the west of the heating maxima. Zonal wind speeds become reslistic when the value of e is 
increased, for example, to 1. For e = 1, the zonal distance of U ma x is just 2° to the west of heating 
maxima. The jet is also more zonally confined and appears symmetric about the peak. Thus, reducing 
zonal wind speed has put a penalty on the location and shape of the jet. The meridional shape is 
independent of e and is thus same as Fig. [3a] Zonal half length and width listed in Table [2] show a 
very short jet for standard and high e values. 


4.2 Gill heating in ideal-physics version of CAM-3.1 

It is generally assumed that the TEJ can be explained with the help of Gill’s formulation. But this 
is quite a simplistic viewpoint since the tropospheric structure has significant diffrences in upper and 
lower halves. More importantly from a dynamical viewpoint, friction is not fully understood and hence 
its parametrization differs significantly in both layers. As was shown before changing friction in G80 
itself changes the zonal wind pattern. 
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Table 2: Maximum heating rates, magnitude and location peak easterly, and jet half-dimension 


Gill model * 


Qmax 


Umax 

(ms 1 ) 


Jet dimension 

6 

( K day 1 ) 

Value 

Lon 

Lat 

Level 

Length 

Width 

0.1 

5.67 

153.48 

11.4°W 

10.5°N 

Top 

34.7° 

15.3°N 

i 


24.15 

2.0°W 


Top 

26.7° 



Gill heating 

in CAM-3.1 ** 




Case 

Qmax 


Umax 

(ms 1 ) 


Jet width 


( K day 1 ) 

Value 

Lon 

Lat 

Press (hPa) 




Qgill_ip_noTeq 


29.87 

65°E 

6°N 

175 

25.1° 


QgillJp_APTeq 

5.67 

25.66 

52.5°E 

6°N 

225 

28.2° 


QgillJp_nG10Teq 


31.16 

72.5°E 

6°N 

175 

19.4° 



* For Gill model length and heating have been dimensionalized by 10° and 4.87 K day 1 


** For Gill heating in CAM-3.1, the 0.5U ma x contour does not close, so jet does not exist 
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The advantage of using an AGCM is that all the necessary non-linearities in the governing dynamical 
equations is present. This along with multiple layers in the vertical will reveal the response to tropical 
forcing which is closer to reality. 




(a) Gill’s solution: U n orm (b) Gill’s solution: heating rate (red contours) and 

Unorm 5 €—0.1, Umax—^ 153.48 m S 



(c) CAM-3.1, ideal-physics: 

Unorm 



(d) Ideal-physics : heating rate (red contours) and 
Unorm 


Figure 3: (a)|(b) Gill’s solution: Unorm contours at Umax zonal location; (c)|(d) ideal-physics simulation 
with analytical temperature profile (APTeq, refer Tabled imposed. Heating rates are in red (contour 
interval 1.5 K day~ x , negative values are dashed), black lines denote U norm contours, ‘cross-diamond’ 
shows the location of maximum easterly wind. Grey region denotes easterlies. 


The same heating profile discussed in section ETTl is imposed in the ideal-physics version of CAM-3.1. 
The peak of zonal heating has has been shifted to 90°E instead of being at the origin as in the Gill 
model. This has been done to facilitate easy comparison with the aqua-planet simulations discussed 
in section IQ Due to zonal symmetry there is be no dynamic change had the source been kept at 
Greenwich meridian. Heating maximizes at 500hPa and becomes zero at the lowest model level and 
tropopause (in CAM-3.1 the tropopause pressure level, TP p , has been defined as TP p = 250—150 cosc/>, 
where <j> is the latitude). Lower tropospheric friction has been retained as before, that is, it is the same 
as HS94. Since upper tropospheric friction has been found to be necessary, the previously described 
Rayleigh friction parametrization has been used. 

Two different temperature profiles have been imposed. The first profile was obtained by zonally 
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and meridionally averaging the 12 th month of AP_90c20n simulation discussed before. Thus a single 
set of points was obtained at each model level. The surface temperature was found to be close to 
298K. The tropopause temperature was fixed as 205K. A second degree polynomial was used to fit 
the values between the lowest model level to the tropopause. Beyond tropopause a linear fit was used 
till topmost model level which was close to 235K. This profile thus has no meridional asymmetry and 
the simulation is named Qgill_ip_APTeq. The second profile used had zonally averaged temperatures 
obtained from July 1 st year noGlOrog simulation. This profile serves to understand the zonal wind 
response in the presence of meridionally asymmetric temperature profiles which actually exist during 
the period when TEJ is observed. This latter simulation is named Qgill_ip_nG10Teq. Additionally 
a set of simulations without any temperatures imposed has also been carried out. This is named as 
QgillJp_noTeq case. This last set is unrealistic and has been done to know the response to heating 
only. The initial conditions are the same as the aqua-planet full-physics simulations. 

The peak easterly zonal wind and its location for all simulations are listed in Table [2] Normalized 
zonal wind profiles at the location of peak zonal wind are shown in Figs. I3cldl The meridional section 
shows that CAM-3.1 qualitatively captures the Gill model features. But the horizontal section shows 
extensive jet lengths. In fact as documented in Table [2J the jet is undefined in the ideal-physics 
simulations. The easterlies extend to the south of the Equator and have greater width compared 
to westerlies. Thus the ideal-physics CAM-3.1 simulation is far from being meridionally symmetric 
despite heating being symmetric. The easterlies also have their maxima between 5°N-10°N. Westerlies 
peak aroundf 20°S. The meridional cross-section reveals that the asymmetry exists at all vertical levels. 
The low-level winds are sometimes four times weaker than those at upper levels. 

In terms of jet magnitude, the value of U max is in close agreement with Gill solution for e = 1 
which represents very high damping in the upper troposphere. The location of U max does not match. 
The zonal distance between U max and heating is between 20°-30°. 

These simulations thus clearly show that the Gill model though very useful is unable to capture the 
true vertical zonal wind structure revealed in a multi-layered model with non-linear interactions and 
more complicated frictional parametrization. It can be expected the when multiple heat sources are 
considered these inadequacies will be further exposed. Further, when full AGCM physics is taken into 
account, there will be further deviations from the Gill model. Beyond crucial and insightful qualitative 
explanations, the analytical formulation of Gill cannot give accurate quantitative answers. 

5 Conclusions 

An AGCM, CAM-3.1, has been run with idealized physics and the simulation has been compared with 
the corresponding full-physics simulations. The ideal-physics simulations used a very simple frictional 


12 



parametrization. In the absence of upper tropospheric friction, using the same heating from full-physics 
simulation caused the jet to be incorrectly simulated. Incorporating upper tropospheric friction made 
the jet closely match the corresponding full-physics simulation. 

The Gill model predicts the presence of a TEJ. But unless very high friction is used, the magnitude 
of zonal winds is very high. Heating from Gill model has been applied in CAM-3.1 to check if an 
AGCM can reproduce a TEJ. The upper tropospheric friction friction which gives realistic simulation 
in CAM-3.1 has been used. It has been demonstrated that Gill heating does not lead to the formation 
of a TEJ. 
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